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Abstract 

In the models with large extra-dimensions, we examine the production of a vector boson (7/ 
Z) in association with the Kaluza-Klein (KK) modes of the graviton via gluon fusion. At the 
leading order, the process takes place through quark-loop box and triangle diagrams and it is 
ultraviolate finite. We report the results for the LHC. We also discuss the issues of anomaly and 
decoupling of heavy quarks in the amplitude. 



1 Introduction 

The problem of large hierarchy between the electroweak and Plank scales can be addressed within 
the models of extra-space dimensions. Among various models of extra-dimensions the model of 
large extra-dimensions, known as the ADD model [1], is one of the very first successful attempts 
in this direction. The predictions of this model can be probed at present day high energy colliders 
Q In the ADD model, the number of space-time dimensions is assumed to be 4 + 5. These 5 
extra-space dimensions are considered flat. The standard model (SM) degrees of freedom live on 
a 3 + 1 dimensional brane, while gravity is allowed to propagate in full 4 + 5 dimensions. To avoid 
any direct conflict with the present day observations, these extra-dimensions are assumed to be 
compact. For simplicity we consider compactification of these extra-dimensions on a torus with 
common scale R/2ir. The fundamental scale of gravity, Ms is related to the 4— dimensional Plank 
scale by 
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1 Though the ADD model solves the hierarchy problem of the SM, it introduces another hierarchy (unexplained 
within the model) between the scale for the size of the extra dimensions (1/R) and the fundamental scale of gravity 
(Ms). But even then it stands as a good phenomenological model to study. 



M% « M 5 S +2 R s . 



(1) 



^From this relation one can argue that if R is large enough the fundamental scale of gravity 
in 4 + 5- dimensions can be as small as few TeV. In 4-dimensions, the 4 + 5 dimensional graviton 
appears as an infinite tower of KK modes with mass spectrum given by, = |^ for the n th 

mode. The interaction of spin-2 components of these KK-modes (h$) with the SM is given by 
an effective Lagrangian 



where T^ v is the energy-momentum tensor for the SM. Although the cross-section for the pro- 
duction of a single KK-graviton is proportional to the inclusive cross-section, obtained by 

summing over all KK-modes, is only l +2 suppressed. Thus if Ms is in the TeV range, it can have 

observable effects at the LHC Since the coupling of each KK-mode with the SM particles is very 
small, the direct production of KK-graviton gives rise to missing energy signal in the detector. 
More details on the model and it's phenomenology can be found in [21 |3l Hj. 

The direct production of KK-graviton (together with QCD NLO correction) has been consid- 
ered by many authors in recent past both at the Tevatron and the LHC [HI El El El EJ- Due to 
large gluon flux at the LHC the gluon initiated processes can be quite important. In this regard, 
we have investigated the KK-graviton production in association with a vector boson (j/Z) via 
gluon fusion. The paper is organized as follows: In section 2 we introduce our process and give 
some details on the structure of the amplitude. Various checks on amplitude and the method of 
computation are described in section 3. Numerical results are presented in section 4. We finally 
conclude in section 5. Some interesting issues related to the amplitude calculation are added in 
the Appendix. 

2 The amplitude 

At the leading order (LO), the process gg —> j/Z + G proceeds via quark-loop diagrams. The 
allowed vertices and Feynman rules for the ADD model are derived in Ref. jl]. Depending upon 
the coupling of KK-graviton with the standard model particles, there are three classes of diagrams. 
The prototype diagram for each class is given in Fig. (1). Other diagrams are obtained just by 
the appropriate permutations of the external legs. There are total of six box and twelve triangle 
diagrams for each quark flavor. However, due to charge-conjugation property of fermion loop 
diagrams, only half of the diagrams are independent. For gg — >■ 7 + G case, we find that the 
amplitudes for diagrams related by charge-conjugation are equal and opposite to each other. This 
implies that at the LO 



This is just an implication of the extension of Furry's theorem in the presence of KK-gravitons, 
photons, and gluons pi)]. By introducing charge conjugation transformation of KK-graviton field 




(2) 



n 



M(g + g ->■ 7 + G) = 0. 



(3) 
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and using charge-conjugation properties of the gluon and photon fields, it can be shown that the 
amplitude does vanish at the LO This result would remain valid to all orders if only QED 
and/or QCD radiative corrections are included. In the presence of the weak interaction, this 
result would not hold to all orders. 

For the case of gg — > Z + G, the amplitude has both vector and axial-vector contributions 
coming from Z-boson coupling to quarks. The vector contribution is similar to 7G case in struc- 
ture. Thus, the above argument when applied in ZG case, implies that at the LO the process 
receives contribution only from the axial part of the Z-boson coupling. For a given quark flavor 
in the loop the amplitude has the following structure, 



raft 



M a p b (g + g^Z + G)=gg 2 s K 



c q A A(m q ,s,t), 



A(m q ,s,t) = A t ri(m g ,s,t) - A box (m q ,s,t). 



(4) 
(5) 



Here g 



cose 'sine J K = V^7rG N ; c\ -- 
and b are color indices for the two gluons. Atm is the sum of the two types of triangle contributions 



|-. s and t are Mandelstam variables. Furthermore, a 
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Figure 1: Prototype Feynman Diagrams for g + g — >■ 7 jZ + G. 



This one-loop process is a LO process and is expected to be free of ultraviolet (UV) and 
infrared (IF) singularities for each quark flavour. We also expect gauge invariance with respect 
to gluon and KK-graviton current couplings. However, due to the presence of anomalies, there 
will be no gauge invariance with respect to axial vector current coupling for individual quark 
flavour. But, since the model is free from anomalies, we expect gauge invariance when we sum 
over all quark flavours. The confirmation of the cancellation of UV and IR singularities and gauge 
invariance with respect to vector and axial-vector current couplings are powerful checks on our 
calculation. We make all these checks as described in the next section. 

2 Thc Graviton field is considered even under charge conjugation as it couples with the energy-momentum tensor 
only. This can be explicitly verified using the ADD model Lagrangian given in Ref. [3]. 
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Of the six quark flavours, we treat the u, d, s, and c quarks as massless. Since the amplitude 
of the process is proportional to T% value, first two generations do not contribute. Therefore, the 
total amplitude, including the contributions from all six quarks, is 

1 fi ab 

E M f(9 + 9^Z + G) = - ^*(— ) [A(m t , s, t) - A(m b , s, t)) . (6) 
<? 

We note that the cross-section is of 0(a^), and therefore this LO contribution can be included 
in, a NLO (pp ^ Z + G + X) mm. 

3 Calculation 

The amplitude for each diagram is written using the Feynman rules given in [I]. However, we 
need not compute all diagrams explicitly. We only need to compute the prototype diagrams. All 
other diagrams can be obtained by suitable permutations of external momenta and polarizations. 
However, due to the presence of 75 in the amplitude, one needs to take extreme care in making 
these permutations in n-dimensions. The permutation should not be across the 75 vertex. Such 
permuted diagrams need to be computed explicitly. Due to the presence of quark-loop, each 
diagram contribution is proportional to a trace over a string of gamma matrices. We compute 
these traces using FORM [TT]. This is the most important part of the calculation. The presence 
of 75 in the trace leads to spurious anomalies in 4-dimensions. More details on this is added 
in appendix. We therefore need an appropriate n-dimensional treatment of 75. We use Larin's 
prescription for 75 to calculate the trace in n-dimensions. According to this prescription 

7m7s = g W^tVt' ( 7 ) 

The fermion loop trace in n-dimensions has an interesting feature. It is independent of the 
fermion mass. After calculating the trace, we express the amplitude in terms of appropriate 
tensor integrals. The box amplitude has rank-four tensor integrals, while triangle amplitude 
has rank- two tensor integrals at the most. At this stage, the amplitude in FORM is converted 
to Fortran routines and is written in terms of tensor integrals and kinematic invariants. The 
Fortran routines to numerically compute these tensor integrals in terms of scalar integrals were 
written in [12] using the methods of Odenbrogh and Vermaseren [13] The reduction is done in 
n = (4 — 2e)-dimensions. The most general form of amplitude after tensor reduction is 

M = J2 ( d ^o) + E ( c ^o) + E MS + E Mo) + K, (8) 

i i % i 

where Ao,Bq,Co and Dq are scalar integrals in the notations of 'tHooft and Veltman [14] . The 
coefficients in front of scalar integrals (also know as scalar form factors) contain information on 
the tensor structure of the amplitude. TZ, known as the rational part of the amplitude, is an 
artifact of UV regularization of tensor integrals. For the massless quark case there is no A term 

3 We have also used an analytical tensor reduction code (in FORM) based on Denner and Dittmaier's work [1 5J , 
to make various checks on our amplitude. 
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in the above expression [^J We verify that the rational term, 1Z being an UV effect is independent 
of the quark mass. 

All the required scalar integrals for the massless quark case are calculated following the method 
of 'tHooft and Veltman [TJ]. We need only the UV and IR singular pieces of these scalars to 
verify the cancellation of UV and IR singularities (which are actually large logs for the case 
of light quarks). We regularize UV divergences (present only in v4 's and -Bo's) in dimensional 
regularization, while IR divergences are regularized by giving small masses to quarks. For the 
case of top quark and bottom quark, we use FF library to calculate the required scalars [16J. The 
above amplitude, after substituting the scalars has following symbolic form 

M = M uv — + M IR 2 ln 2 (m 2 ) + M IR ln(m 2 ) + M F . (9) 

Mp is the finite part of the amplitude and includes the rational term, 1Z. Due to very large and 
complicated expression for the amplitude, we compute the amplitude numerically before squaring 
it. This requires computation of polarization vectors for gauge bosons and for the graviton. We 
have chosen helicity basis for them. It also helps in making additional checks on our calculation by 
verifying relations among helicity amplitudes. The polarization tensor for graviton is constructed 
from polarization vectors of massive vector bosons as suggested in [3]. 

As discussed in the previous section, the process is expected to be finite. We verify these 
expectations for our amplitude. Both the massive and massless contributions are UV finite. We 
observe that each triangle diagram is UV finite by itself, while Box amplitude is UV finite after 
adding all the box contributions. Fermion loop diagrams are known to be IR finite in both 
massive and massless fermion cases for any kind and any number of external particles attached 
to the loop [T7]. In the massless quark case, we check that each diagram is IR finite and therefore 
M/#2 = Mir = 0, holds for the full amplitude. In the case of triangle class of diagrams involving 
graviton-gauge boson coupling, we choose £ = 1 (the Feynman gauge) for gluons case and £ = oo 
(the Unitary gauge) for Z-boson case. As expected the calculation does not depend on any specific 
choice of the gauge parameter. Finally, we check gauge invariance of the amplitude with respect 
to the two gluons by replacing their polarizations with their respective momenta. We observe that 
(only after using 75 prescription in the trace) both the massless and massive quark contributions 
are separately gauge invariant with respect to the two gluons. We also check gauge invariance 
with respect to the Z-boson. Because of the anomaly, the two contributions are not separately 
gauge invariant but the total amplitude is gauge invariant up to the t-quark mass as expected 
due to the explicit breaking of the chiral symmetry. All these checks on our amplitude have been 
made both numerically as well as analytically. The issue of anomaly in the amplitude is discussed 
in appendix. 

4 even for the massive case, the Aq can be replaced by, Ao(m 2 ) = m 2 (l + Bq(0, m 2 , m 2 )). 
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4 Numerical Results 



In this section, we present the results of our calculation for the process pp — > Z + G + X, at 
the LHC. The results would depend on the two parameters of the ADD model - (i) the number 
of extra-space dimensions, 8 and (ii) the fundamental scale of gravity, Mg. We explore this 
dependence and other features of the process. In Fig. (1), we have plotted the cross-section as a 
function of the center of mass energy. Here we have chosen 5 = 2 and M s = 2 TeV. In addition, 
we have applied following kinematic cuts: 

P# > 30 GeV, \r] z \ < 2.5, < M s . 

The cut on partonic CM energy or equivalently on ZG invariant mass is known as truncated 
scheme in literature. This cut is related to the fact that theoretical predictions within the ADD 
model, which is an effective theory, can be valid only below the fundamental scale, Ms- It will be 
interesting to probe the sensitivity of our predictions on this kind of constraint. We will further 
comment on the issue at the end of this section. We have also chosen factorization and renor- 
malization scales as /// = /j,r — E% (= a/ Mf + (P/) 2 ). We have used NLO parton distribution 
functions CTEQ6M. We note that at typical LHC energy, the cross-section would be of the order 
of few fb. These cross-sections are much smaller than expected. The smallness of the cross- 
section is due to two-orders of magnitude cancellation in the amplitude between the box-type 
and triangle-type diagrams. This destructive interference occurs due to the relative minus sign 
between the box-type and triangle-type diagrams contributions. Still one may expect few hundred 
of such events after a few years of LHC operation at 14 TeV CM energy. 

We now illustrate the various aspects of this process at the center of mass energy of 14 TeV. 
The parton distribution functions, factorization/renormalization scales, and the kinematic cuts 
are as above. 




V(S) (TeV) 

Figure 2: Cross-section as a function of collider energy at the LHC. 

In Figs. (3) and (4), we show the dependence of total cross-section on the two parameters 
of the ADD model, M s and 5. As 5 or M s is increased, the density of states for KK-modes 
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Figure 3: Dependence of total cross-section on number of extra dimensions 5, for M s = 2 TeV. 




Figure 4: Dependence of total cross-section on the scale M s , for 8 = 2 
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Figure 5: Transverse momentum distribution of Z for M s = 2 TeV and 5 = 2. 



Figure 6: Rapidity distribution of Z for Ms = 2 TeV and 5 = 2. 
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Figure 7: KK-Graviton mass distribution for Ms = 2 TeV and 5 = 2. 



falls and therefore the cross-section also goes down. The kinematic distributions of the transverse 
momentum and rapidity of Z are plotted in Figs. (5) and (6), respectively. Fig. (7) shows the 
KK-graviton mass distribution. Clearly, the cross-section is peaked at a relatively larger value of 
about 200 GeV. We have also calculated the total cross-section for a cut of 400 GeV to avoid 
the SM background as suggested in [2J. We find the cross-section of about 0.2 fb and it's almost 
10% of the NLO QCD correction calculated in 018]. 

Next, we study the scale dependence of the total cross-section. We vary the factorization 
scale around it's central value, \ip = E%- The cross-section varies by about 30% by changing \ip 
in the range between E^/2 and 2E^. The uncertainty in our calculation due to the choice of 
PDFs is between 5 — 20%. This is estimated by computing cross-section for other CTEQ parton 
distributions - CTEQ6D, CTEQ6L1, and CTEQ6L2. We have also studied the effect of taking 
large fermion mass limit in the loop. For a given phase-space point, we vary the top-quark mass 
and observe that the amplitude-squared (which includes both the bottom and top quark contri- 
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butions) approaches a constant value as m t — > oo. This implies the complete decoupling of the 
top quark, i.e. the top quark contribution of the amplitude goes to zero in large m t limit. It is 
expected from the decoupling theorem [18J. This feature has been plotted in Fig. (8). The change 
in slope around m t = ^ corresponds to physical threshold after which top-quark pairs can't be 
produced on-shell and the amplitude is real. 

16 
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10 

% 8 
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4 
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m, (GeV) 

Figure 8: Decoupling of t-quark in m t — > oo limit. 

Finally, we comment on the UV sensitivity of our theoretical results. As mentioned before, we 
have presented all the results for which the partonic energy is below the fundamental scale 
Ms, i.e. the truncated scheme. It has been argued in Ref. [2] that, if relaxing this constraint 
does not change the results significantly, then results of the effective field theory can be trusted. 
We find that the total cross-section in two cases (with and without cut on differs by 20% 
for 5 = 2 and Ms = 2 TeV. This difference, as expected, increases for larger values of 5, while 
decreases with increasing Ms- 

5 Conclusion 

We have presented the leading order calculation for the production of a KK-graviton in association 
with a vector boson via gluon fusion in the ADD model. Results have been presented for cross- 
sections and some important kinematic distributions at the LHC We find that the amplitude for 
the jG case is zero at the LO, while in ZG case only axial part of the amplitude contributes. We 
have studied the qualitative features of our amplitude with respect to the anomaly structure and 
the decoupling of heavy fermions in the loop. We saw that both the triangle and box diagrams 
give rise to anomaly in the amplitude and this anomaly does cancel after considering contributions 
of other quark flavors. We also observed that the amplitude for a given quark flavor in the loop 
vanishes in the large quark mass limit verifying the decoupling theorem. 

The total cross-section for the gluonic contribution is small compared to the NLO results for 
the ZG production at the LHC reported in [TJ [S] . The smallness of the cross-section is due to 
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the destructive interference between the box-type and triangle-type diagrams. There is a relative 
negative sign between the two contributions. This leads to the two-orders of magnitude cancella- 
tion that takes place between box and triangle contributions at the amplitude level. Still, there 
can be about 100 events per year when LHC will operate with peak CM energy and luminosity. 
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6 Appendix 

6.1 Issue of 75 and Chiral Anomaly 

The triangle amplitudes for our process have VV A-structuie, while box amplitude has VVVA- 
structure. The issue of chiral anomaly is well known in linearly divergent fermion loop triangle 
diagrams with VIM-structure. Although box diagrams have WIM-structure, they are also 
linearly divergent because of the nature of graviton coupling with quarks. Therefore, both the 
triangle and box diagrams will give rise to anomalous contribution. The presence of anomaly 
affects the gauge invariance of the amplitude. For reliable predictions, our amplitude should be 
gauge invariant with respect to all the currents. Let's try to understand the issue of anomaly at 
the level of general structure of one loop amplitudes given in eqn. (8). 

As discussed earlier, the rational term TZ in eqn. (8) is due to the UV regularization of ten- 
sor integrals and it is independent of masses of internal lines. So, in dimensional regularization, 
whenever \ UV-pole hits an 0(e) term, we get a finite contribution. TZ is collection of all such 
terms obtained during the reduction of tensor integrals in n = 4 — 2e dimensions. It should be 
clear from the general structure of one-loop amplitudes, that in a gauge invariant amplitude TZ 
should be separately gauge invariant. In the case of one- loop calculation of fermion- loop diagrams 
in dimensional regularization, the issue of anomaly is related to the presence of 75 in a linearly 
(or worse) divergent integrals. There is always an ambiguity in carrying out full calculation in 
dimensional regularization, as 75 is strictly a 4-dimensional object. But this ambiguity arises 
only in fixing the rational part of the amplitude. Any 0(e) structure of 75, if it at all exists in 
n = 4 — 2e dimensions, will contribute to the rational term as explained above. Thus anomalies 
affect only the rational part of the amplitude. They are also independent of the fermion mass. 
In the presence of anomalies, the gauge invariance of the amplitude is violated by these rational 
terms only. 

The anomaly being associated with 75 should spoil only the axial-vector current conservation 
and the conservation of vector current must hold. But it is not very clear if this feature can 
always be ensured by treating 75 in 4-dimensions. The 4-dimensional 75 may lead to spurious 
anomalies in the amplitude, resulting in non-conservation of vector currents. To ensure the vector 
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current conservation in the amplitude, it is advisable to use an appropriate prescription for 75 
in n-dimensions. It is important to note that if we don't treat 75 appropriately in n-dimensions, 
even the relation between charge conjugated fermion loop diagrams does not hold and as expected 
the violation is only in the rational terms. Thus by using an appropriate prescription for 75 in 
n-dimensions we generate a rational term consistent with various symmetries of the amplitude. 
The 75 prescription given in eqn. (7) does this job. One should be careful using this prescription 
when generating other amplitudes by permuting external momenta and polarizations in prototype 
amplitudes. Any amplitude in which position of 75 should also change under such permutations 
should be computed separately. Thus within this prescription we require two prototype box and 
four prototype triangle (two for each class) amplitudes to generate the full amplitude. 

The amplitude thus obtained for a given quark flavor in the loop is anomalous with respect 
to the Z-boson current only. As we said earlier this anomaly in axial-vector current should also 
vanish. Since chiral anomaly is independent of quark mass and our amplitude is proportional to 
the T 3 value of each quark flavor, the anomaly in axial-vector current does vanish when contribu- 
tions from all six quark flavor is considered as in eqn. (6). 

In the light of above understanding of anomalies, it is clear that even if we don't use any pre- 
scription for 75, as far as total amplitude in eqn. (6) is concerned, we must get a gauge invariant 
result with respect to all the currents. Although amplitude for each quark flavor will not be gauge 
invariant any more, their difference in eqn. (6) is always gauge invariant. This is because the 
rational terms being independent of quark mass cancel. We have verified this in a separate calcu- 
lation and all the numerical results described in section 4 agree with this way of doing calculation. 



6.2 Decoupling of heavy fermions and Rational term 

In section 4, we found that if we use n-dimensional prescription for 75, the amplitude for a given 
quark flavor in the loop goes to zero as the large quark mass limit is taken. The decoupling 
of heavy fermions is, in general, expected in any UV finite fermion loop amplitude involving 
couplings which are not proportional to the fermion mass in the loop. This is a very nice feature 
of fermion loop amplitudes. It can be utilized to fix the rational term of the amplitude without any 
ambiguity. It is particularly important in our calculation where the generation of correct rational 
term requires an n-dimensional prescription for 75 consistent with the vector current conservation. 
Since our amplitude is UV finite and decoupling of heavy fermions does hold, it implies that there 
must be a cancellation between the mass independent rational term and non-rational part of the 
amplitude in eqn. (8), as large fermion mass limit is taken. Therefore, we do not need to calculate 
rational terms explicitly. The rational term 1Z is simply the negative of the non-ration part of the 
amplitude in large fermion mass limit. Though decoupling does not hold if Higgs is attached to 
the fermion loop, the structure of the fermion loop amplitude (with Higgs) suggests that one can 
still fix the rational term by taking large fermion mass limit. For that one should take away the 
overall (explicit) factor of the fermion mass before taking the limit. Once again the amplitude 
goes to zero indicating cancellation between rational term and non-rational part of the amplitude. 
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